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1. Introduction 

Recent years have witnessed significant improvements in charmed-meson leptonic and semilep- 
tonic decays, both in experimental measurements and in calculations of the hadronic transition am- 
plitudes with lattice QCD. A puzzle has arisen, namely a discrepancy of approximately 3.5a in the 
rate of the leptonic decay D s — > Iv, where I is a muon or % lepton If the measured counts have 
not fluctuated high, and the lattice QCD calculations are confirmed (by further calculations with 
2+1 flavors of sea quarks), then this may be a signal of physics beyond the Standard Model [gj. 

If non-Standard interactions mediate cs — ► vl, then they also alter, at some level, the rate and 
^-distribution of D — > K/J.V. (D — ► Krv is kinematically forbidden.) In this paper, section || recalls 
the origin of the leptonic discrepancy, incorporating new, preliminary results. Section [| updates the 
new-physics analysis of Ref. ^ and extends it to encompass semileptonic decays. Then section |] 
discusses the phenomenology of semileptonic decays in the context of new physics. For lattice 
QCD the main conclusion, discussed in section ||, is that precise calculations of the semileptonic 
form factors, including a tensor form factor defined below, are vital. 



2. Leptonic Decays 

In the Standard Model the partial width for D s —> £ve is 

T(D S - lv t ) = n ^fl s \ GpV* s m(\ 2 (1 - mj/m 2 D f, (2.1) 

where the decay constant fo s is defined by (0\sy^Y5 c \^s(p)} = ifD s P^, and is also computed via 
(m c +m s )(0\sY5c\D s (p)) = —ifcm 2 Ds ; PCAC ensures that the two definitions are the same. The par- 
tial widths are small: for muonic decays owing to the helicity-suppression factor m^; for T-leptonic 
decays owing to the phase-space factor (1 — m 2 /mj)) 2 . Experiments measure the branching frac- 
tion B = Ttd s but usually quote fo s assuming that no non-Standard amplitude contributes to T. 

In this sense, fu s has been measured recently by the BaBar ||3|], Belle [Hp, and CLEO Col- 
laborations. The experiments measure B(D S — > Iv) directly, without complicated modeling of the 
events or background, and the experimental errors are principally statistical. Radiative corrections 
are at most 1-2%, and the discrepancy cannot be explained with any value of |V„ | consistent with 
a unitary n x n CKM matrix In summary, it seems sound to take the experimental measure- 
ments of fo s at face value, yielding Table |]. Treating both statistical and systematic uncertainties 
in quadrature, the average of the measurements in Table || is 

/ D Jexptavg = 272±8MeV (2.2) 

combining jUV and TV and including new results reported by CLEO at conferences through Septem- 
ber 2008 [01. Separate averages for the two final states are in Table |. 

Now let us turn to lattice QCD calculations of fj) s . There are two calculations with 2+1 flavors 
of sea quarks, the first from the Fermilab Lattice and MILC Collaborations [Jsj] and more recently 
from the HPQCD Collaboration [g]. These are 

fD. Ihpqcd = 241 ± 3 MeV, f Ds | Fermilab . MILC = 249 ± 1 1 MeV, (2.3) 



where the Fermilab-MILC result is an update presented at this conference by Mackenzie [10]. Both 



calculations use the improved staggered Asqtad action for the sea quarks, taking advantage of the 
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Table 1: Recent experimental values of fo s . The preliminary update from CLEO can be found in Ref. [Q. 



final *itntp 






(MeV) 






end 2007 


2008 update 




BaBar [j| 


283 ±17 ±16 




fiv 


CLEO [@1 


264±15± 7 


265.4±11.9±4.4 


fiv 


Belle [@] 


275 ±16 ±12 




TV (T — > 7TV) 


CLEO [§ 


310±25± 8 


271±20±4 


tv (t — > evv) 


CLEO [§ 


273±16± 8 






our average 


273 ±11 


271 ±10 


TV 


our average 


285 ±15 


272 ±13 



freely available MILC ensembles [11]. In the following, I use the average of the two results in 
Eq. ( p3| ) with no correlation, because the dominant systematic errors differ. 

The experimental average of fo s lies 12.5% above that of lattice QCD, and the significance of 
the discrepancy is 

3.5a = 2.9a ©2.2a, (2.4) 

where the two entries on the right-hand side are for ptv and TV separately. Before the update 
from CLEO [0] the discrepancy was 3.8a = 2.7a ©2.9a [g]. Reference [jl|] omits BaBar's result 
and reports 3.4a, and with CLEO's update this approach yields 3.2a. One should bear in mind 
that the yardstick for a is the experimental statistical error: Were one to double HPQCD's error 
(without justification), the total discrepancy would remain 3.0a. Indeed, the same methods agree 
with experiment for f % , fx, charmonium mass splittings, and especially mo s , m D +, and f D + [^, O]. 



3. New Physics 

If the discrepancy cannot be traced to a fluctuation or error in either the measurements or the 
calculation(s), then one should turn to non-Standard physics as an explanation. The new particles 
must be heavy to have escaped direct detection, so one may consider an effective Lagrangian 

ifeff = M- 2 Ci(sf r5 c)(v eL r^ L ) +M- 2 C e P (s- 75 c)(v iL £ R ) -M- 2 C e v (sfc)(v iL y^ L ) 

+ M- 2 Cl(sc)(v fL e R ) +M- 2 C e T (so^ v c)(v eL o tlv l R )+H.c., (3.1) 

where M is a high mass scale. This Jzf e ff extends the effective Lagrangian of Ref. ^ to include 
interactions that mediate D — ► Kiv. The experiments do not identify the neutrino flavor or helicity, 



but Eq. (3.1) assumes v^ L to be of lepton flavor I and omits right-handed neutrinos. In this way 
the resulting non-Standard amplitudes can interfere with the Standard IV-mediated amplitude and 
explain the sought-after effect of 10-15% (in the amplitude). 



These effective interactions change the rate for leptonic decay, by substituting in Eq. (2.1) 



G F V*m ( -> G F V*mt + -p^mf + —=£■ ^i— . (3.2) 

cs cs V2M 2 V2M 2 m c + m s 

Because (conventionally) V cs is real, one sees that one or both of Cr A , C l P must have a positive real 
part. If only one of these reduces the discrepancy to la, one can derive the bounds 
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M{ReC e A y l/2 < 855 GeV, M(ReC e p y l/2 < 1070 GeVy/m T /m e , (3.3) 

updating Ref. [§] to reflect CLEO's new preliminary measurements and treating the jj,v and TV 
discrepancies as a single effect. 

The effective Lagrangian can arise from the tree-level exchange of non-Standard particles, in 
which case M is simply the new particle's mass. Reference [^] found a few possibilities. One 
is the ^-channel annihilation through a charged Higgs boson, in a new model designed so that the 
Yukawa couplings satisfy y s <C y c and y c , y r ~ 1 . But this model also has y c j < y s , thereby predicting 
a 10-15% deviation in the amplitude for D + —* £ + v. This is now disfavored, because CLEO's new 
measurement of f D + [13] agrees perfectly with lattice QCD ^ |9[ [Kj. Another candidate is the 



^-channel exchange of a charge +| leptoquark, which can arise in various ways, all of which 
are disfavored by non-observation of T — ► }iss. The most promising mechanism is the M-channel 
exchange of an SU(2)-singlet, charge — A leptoquark, namely a particle with the quantum numbers 
as a down-type scalar quark d in supersymmetric models. It couples via the /^-violating Lagrangian 

-2lq = *u {cdl - s L v c a ) d+ 4e c R f R d + H.c. , (3.4) 

where the superscript c denotes charge conjugation, and fc^ and are complex parameters (in 
general, entries of Yukawa matrices). When M = mj » mo s one can derive «Sf e ff with 

C A = c v = il^l 2 , Cp = Cg = l^2t^2£ = — 2Cf . (3.5) 

If K2e is independent of £ and either x mi or \t^/tcit\ <C mim c /mj), then these interactions 
could explain why the discrepancy appears in both jitv and TV channels. Generalizations of 
Eq. (3.4) appear in non-Standard models that modify the interference phase of B®-B® ex- 
plain quark masses Jl5|], induce deviations in — > Iv Jl6|], generate neutrino masses [|l7|], or 



enhance rare D decays [18] 



4. Semileptonic Decays 

To obtain further information about a possible non-Standard cause of the effective scVl vertex, 
one can turn to other processes. One would be the production charmed quarks in neutrino scattering 
off strange sea quarks in nucleons. Another set consists of the semileptonic decays D — > K pL + v^, 
D + — ► K }X + v^, and their charge conjugates. A full understanding of these decays will require 
lattice QCD calculations of the hadronic transition. 

Let us start by reviewing the kinematics of three-body decays. Let the D-meson, kaon, lepton, 
and neutrino 4-momenta be denoted p, k, i, and v. There are two Lorentz independent invariants, 
which may be taken to be E( = p ■ £/mjj and Eg = p ■ k/mj), namely the lepton and kaon energies 
in the D meson's rest frame. Often instead of Ek the mass-squared of the leptonic system is used, 
q 2 = m 2 D + m\ — ItudEk, q = £ + v= p — k. For brevity the formulae given below use both Eg 
and q 2 . The kinematically allowed region is shown in the Dalitz plot, Fig. [T[ The discussion given 
below is somewhat simpler with the variable 

E^ = P —- P -*t± = Ei -l(m D -E K )(l + mj/q 2 ) , (4. 1) 
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and the allowed region, for fixed E K is -£™ ax < E ix < E™ x , Ej%* = ±(1 - mj / q 2 ) J E 2 - m\. 
The allowed region for Ek is mg <E K < (ml + m 2 K — mj) /2mjj, or mj < q 2 < (mo — m^) 2 . 
The doubly-differential rate for D — > K£v is 



+ 



+ 



2 2 

T — rrq 
Am 2 D 



m e [G F V*+G\ 



m D -m K 



+ G 



t m D -m K 



m r — m.. 



\h{q 2 



Am 2 D 



(E 2 -m 2 K )(\-mj/q 2 ) + ^-E i 



\G e T \ 2 \f2(q 2 )\ 2 



(4.2) 



^ (E 2 -ml) (I- mj/q 2 ) Re [{g f V* + G* v ) Gpf + (q 2 )f^ (q 2 ) 
2m$ 



E i± Re 



m D 



2q 2 
mi, 



mAG F V* + G\ 



ntn — m 



D " l K i f~»t D 



mn-mt 



m r — m K 



G F W cs ^G^)h(q 2 )fl(q 2 ) 



m £ [G F V*+Gi 



m D -m K 



+ G 



e m D- m F 



m r — m<. 



G l T *Mq 2 m(q 2 ' 



where G VST = C vs T /V2M 2 , and the form factors / + , /o, and /a are defined via 

(K(k)\sfc\D(p)) = ^+^-^^q^f + (q 2 ) + ^^q^Uq 2 ) 
(K(k)\so» v c\D(p)) = im D l ( P n v -p v k»)f 2 (q 2 ), 



(K(k)\sc\D(p)) 



m D~ m K 



M 



(4.3) 
(4.4) 
(4.5) 



m c - m s 

and fo appears for both the vector and scalar currents owing to CVC. Integrating over lepton energy 



dT 

dEv 



+ 



(2w)3 

±_ 

Am 2 D 



2 2nV+'« ; 



3</ 2 



i\G F V; s + Gi\ 2 \f + (q 2 )\ 2 



m ( [G F V c * s + G v 



m D -m K 



+ G 



m r — ra, 



\fo(q 2 



(4.6) 



+ (El-ml)^^- \G £ T \f 2 (q 2 )\ 
5m D 

-2^(E 2 - m\) Re [ (g f V* + ) Gpf + (q 2 )f 2 (q 2 ) 
mjj L V / 

Note that two terms with interference between form factors vanish after integration: the vari- 
able E(± renders this feature especially transparent. We shall use these formulae to diagnose how 
new interactions mediating D s — ► £v would alter the semileptonic rate and differential distributions. 

In many non-Standard models, including those with the charge — j leptoquark d, one finds 
Cy = C A and C| = Cp (cf. Eq. (|3~3|)). If the fo s puzzle is solved by the C l A interaction, then the C v 
interaction generates a similarly large enhancement in the semileptonic rate. This consequence is 
easily seen from the \f+\ 2 contribution to the rate: the other Standard contribution, with \fo\ 2 , is 
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suppressed by (m^/m^) 2 , which is 3 x 10~ 3 (7 x 10~ 8 ) for \i (e). On the other hand, if the fo s puz- 
zle is solved by the Cp interaction, then it will be difficult to observe the companion semileptonic 
contribution. To enhance both D s — ► TV and D s . — > /iv, some mechanism should lead to Cp oc mi, 
and then (in the leptoquark example) C| °< mi, Cj °< also. In that case, the non-Standard contri- 
butions are small corrections to a suppressed contribution. 

The doubly-differential rate suggests a challenging way to observe the effects of non- vanishing 
C s and Cf. In the asymmetry 

_ r(E e± > 0) - r(E e± < 0) _ N(E t± > 0) -N(E e± < 0) 
1 T(E l± > 0) + r(£ £± < 0) N(E e± > 0) < 0) 

everything but the last two lines of Eq. ( |4.2| ) cancels. To obtain a 7% measurement of gf±, one 
would need around 10 7 semimuonic events in each half of the modified Dalitz plane. One can 
generalize this asymmetry to any region in the E K -E{± plane that is symmetric about E#x = 0, or to 
any moment of the distribution odd in E&±. For example, when the kaon momentum is low, phase 
space naturally suppresses the |/+| 2 contribution, perhaps helpfully. 

5. Conclusions 



From Eqs. ( |4.2| ) and ( |4.6D one sees that the first concern of future lattice calculations is to 
improve on the 7% uncertainty of the only 2+1 flavor calculation of f+(q 2 ) [19]. With current 



semielectronic measurements [20, 21], one could test for new contributions to the scv e e vertex (for 



which there is not yet any evidence). Semimuonic measurements are needed for a direct test of the 
fo s puzzle. Once the event yields become high enough to measure sfy_, it will be necessary to have 
accurate calculations of the scalar and tensor form factors, /o and fa. CLEO and the B factories 



have somewhat more data to analyze, and BES-III should record thousands of events [22]. 
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Figure 1: Dalitz plots for D — > K£v, with mg = m,u. The left panel shows Ef, the lepton energy in the D 
meson's rest frame; the right panel shows E(±, defined in Eq. (4.1). The dashed (blue) lines show E(± = 0. 
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A. All Semileptonic Formulas 

To treat the SM and NP efficiently, we shall now write the effective Lagrangian as 



where 



Jz? eff = M- 2 Ci(sf75c)(v L r^ L ) -M- 2 4(^c)(v l7m 4) +M- 2 C e P (sy 5 c)(v L e R 
+ M- 2 C l s (sc)(v L £ R ) +M- 2 C i T {sa^c){v L a llv t R ) +H.c, 

c" VA = V2M 2 G F v; s +ci A , = g f v; s + 



(A.l) 



(A.2) 



V2M 2 cs V2M 2 ' 

This Jz? e ff mediates D — > Klv. (By analogy, one can extend this to the semileptonic decay of any 
pseudoscalar meson.) Let the D-meson, kaon, lepton, and neutrino 4-momenta be denoted p, k, £, 
and v, as above. The amplitude is 



(£vK\i^ eS \D) = iM- 2 u(v)\{\ + y 5 ) C e s {K\sc\D) - C^(K\sy^c\D 



(A3) 



+ C T a^ v (K\sa^ v c\D) v(l). 



The hadronic matrix elements are re-expressed as form factors in Eqs. d4.3|)-([4~4l). For the leptons 
we require the spinor combinations 

u(v)±(l + y 5 )(t+fi x v(£) = 2fi(v)^(l-fc)v(*), (A.4) 
"(v) 5(1 + 75)^) = -mtu(v)l(l + y 5 )v(£), (A.5) 
(p"k v -p v k^u(v)^(l + y 5 )io^v(£) = -2u(v)tJ$(l + K)v(£), (A.6) 



where for any four-vector r, = — (r-qjq 2 )q^. Inserting Eqs. (4.3)— (4.4) and (A..4)-(A.6) into 
Eq. (B, 



(£vK\i^ eS \D) = iM~ 2 <j ( mtOf - 



mf, — m 



K , n l m D~ m K 



m c — m s 

,2\ r.M-\ 



w(v)I(l + 7 5 )v(£)/o(<? 2 ) (A.7) 



2C^(v)^i(l-y 5 )v(^)/ + (^)-2C^ 1 «(v)^i(l + r 5 )v(^)/ 2 (^ 



The differential rate is (see PDG) 

d 2 F 1 1 



dE^dEi (2k) 3 Smo 



£ \(£vK\iJ? ef[ \D)\ : 



(A.8) 



spins 



where Ek = p -k/mo and E( = p -i/mo are the energies of the kaon and lepton in the rest frame of 
the D meson. To sum over lepton and neutrino polarization states, we need 



£ «(v)i(l + y 5 )v(£)v(£)i(l - YsHv) = 2v-£ = q 2 



III, ■ 



spins 



£ M -(v)^I(l-7 5 )v(W)i(l + 75)^«(v) = -p 2 Aq 2 -m 2 )-(2p ± 



(A.9) 
(A. 10) 



spins 



£ M(v)^i(l + 7 5 )v(£)v(£)i(l-75)# ± «(v) = f{2p ± .e) 2 -n$pitf-n$), (A n) 



spins 
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£ w(v)±(l + 7 5 )v(W)i(l + y 5 )/^(v) = 2mtp ± -e, 

spins 

£ u{v)\{\ + y 5 )v{l)v{i)\{\-y 5 )^ L u{v) = -2q 2 P± -£, 



spins 



spins 



(A. 12) 
(A. 13) 
(A. 14) 



in which p\ = — {E\ — m\)ml/q 2 . Note that E\ — m\ is nothing but the kaon's squared three- 
momentum in the D-meson rest frame. It is convenient to change from the variable Eg to E£± = 
p± • £/mjj, defined in the text. With this variable, it is easy to see which contributions vanish after 
integrating over lepton energy. 

The doubly-differential rate is then 



d 2 T 



1 m D 



dE K dE l± {2nf 2M 4 



+ 



+ 



2 2 

q —m\ 
Am 2 D 



(E 2 K -m 2 K 



- t ml - m\ 



2 q -m e 



-4E 



q L 



2\|2 



K \f + (q Z )\ 



meC' v y 



+ C 



7 m D — m K 



m c — m s 



2\i2 



1/0(^)1 



m h F 2 n?\ q ^ I F 2 
Y \ K m K) ~2 r— t L 

m D q z 

Efj^Re 

m D 



2\|2 



\M<t)\ 



Re 



2\ t*l„2^ 



2 2 2 2 

m fl m D~ m K , n i m D~ m K \ f,t* f („2\ f * („2\ 



m r — m K 



+ 



2i 



ni(C v 



mn-mi 



m c — m s 



(A. 15) 



The singly-differential rate is 



dT 



1 mo 



dE K (2tt) 3 2M A 



ml 1 



+ 



Ami 



mfO 



rrin — m v 



mt 



, m^, — m 



(E K - m K 



2 ,2q 2 + m\ 



3q 2 



141 \U{q 



2\|2 



D " l K 



m c — m s 



2\|2 



\fotf)\ 



?2 „.2W 2 + 2ffl 



+ {Ek-O- , 2 



(A. 16) 



(■ vW)i 2 - ^(4-^l)Re [447 + (<7 2 )/ 2 V) 



= ^y^i(i-f) 2 {(4-».i.) 



2 2 2q 2 +m 2 



3q 2 



G F V* + 



V2M 2 



2\|2 



\M<f)\ 



+ 



Am 2 D 



me(G F V* + 



a 



V2M< 



+ (4-1)^ 



3mf 



_ 2 ^i(4-m 2 )Re 



V2M 2 



m 2 D -m 2 K 



\f2(q 2 )\ 2 



C| m 2 -m| 



\/2M 2 m c -m s 



2\\2 



1/0(^)1 



(A. 17) 



V2M 2 / V2M 2 



9 



